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Abstract 

We report results of searches for charmless hadronic B meson decays to pseudoscalar (tt^, 
K ± , n° or K$)-vector(p, K* or uS) final states. Using 9.7 x 10 6 BB pairs collected with 
the CLEO detector, we report first observation of B — ► n~p°, B° — ► 7r ± p T and B~ — > 
which are expected to be dominated by hadronic b —>■ u transitions. The measured 
branching fractions are (10.4±|;1 ± 2.1) x 10" 6 , (27.6±f;| d= 4.2) x 10~ 6 and (11.3±|;| ± 1.4) x 
10~ 6 , respectively. Branching fraction upper limits are set for all the other decay modes 
investigated. 

PACS numbers: 13.20.He,13.25.-k,13.25.Hw,13.30.Eg,14.40.Nd 
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CP violation in the Standard Model (SM) is a consequence of the complex phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix [M. The study of charmless hadronic decays of 
B mesons plays a key role in testing the SM picture of CP violation. For example, the angle a 
= arg [(— V^V^) / (y u dV* b )] of the unitarity triangle can be measured by performing a full Dalitz 
analysis of the decays B°(B°) — > 7r + p~, tt~ p + and 7r°p° 0. While the CLEO data do not yet have 
the sensitivity for the CP violation measurements, experimental measurements of these decay 
modes will be useful to test various theoretical predictions that typically make use of effective 
Hamiltonians, often with factorization assumptions (3|. Recently, it has been suggested 0, with 
model dependency, that published experimental results on charmless hadronic B decays indicate 
that cos 7 < 0, in disagreement with current fits to the information most sensitive to CKM matrix 
elements [||. 

In this Letter, we present results of searches for B meson decays to exclusive pseudoscalar- vector 
(B — > PV) final states that include a pseudoscalar meson 7r , 7T° or K® and a vector meson 
p, K* or to. In particular we present first observation of the decays B~~ — > n~p°, B° — > 7r ± / o =F and 
B~ — > it~lj (charge-conjugate modes are implied) which are expected to be dominated by hadronic 
b — > u transitions. Our results supersede previous CLEO results on these decay modes 

The data were collected with two configurations (CLEO II g and CLEO II. V |J) of the CLEO 
detector at the Cornell Electron Storage Ring (CESR). They consist of 9.1 fb _1 taken at the 
T(4S), which corresponds to 9.7 x 10 6 BB pairs, and 4.4 fb _1 taken below BB threshold, used 
for continuum background studies. The data sample contains a factor of 3 more statistics than 
previously published results ||. In addition, the CLEO II data were reanalyzed with improved 
calibration and track-fitting, allowing for larger geometric acceptance and more efficient track 
quality requirements. 

The final states of the decays under study are reconstructed by combining detected photons 
and charged pions and kaons. The detector elements most important for the results presented here 
are the tracking system, which consists of several concentric detectors operating inside a 1.5 T 
superconducting solenoid, and the high-resolution electromagnetic calorimeter, consisting of 7800 
CsI(Tl) crystals. For CLEO II, the tracking system consists of a 6-layer straw tube chamber, a 
10-layer precision drift chamber, and a 51-layer main drift chamber. The main drift chamber also 
provides a measurement of the specific ionization loss, dE/dx, used for particle identification. For 
CLEO II.V the straw tube chamber was replaced by a 3-layer, double-sided silicon vertex detector, 
and the gas in the main drift chamber was changed from an argon-ethane to a helium-propane 
mixture. 

The resonances in the final state are identified via the decay modes p — ► tttt, K* — ► Kir 
(K*° — ► K + tt~, K* + — ► K + tt°) and u — > 7r + 7r _ 7r°. Reconstructed charged tracks are required to 
pass quality cuts based on their track fit residuals and impact parameter in both the r—<f> and r~z 
planes, and on the number of main drift chamber measurements. Each event must have a total 
of at least four such charged tracks. The dE/dx measured by the main drift chamber is used to 
distinguish kaons from pions. Electrons are rejected based on dE/dx information and the ratio of 
the measured track momentum and the associated shower energy in the calorimeter. Muons are 
rejected by requiring that charged tracks penetrate fewer than seven interaction lengths of material. 
Pairs of charged tracks used to reconstruct Kg (via K s — > tt + tt~) are required to have a common 
vertex displaced from the primary interaction point. The invariant mass of the two charged pions 
is required to be within two standard deviations (a) of the known K s mass [1C]. Furthermore, the 
K s momentum vector, obtained from a kinematic fit of the charged pions' momenta, is required to 
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point back to the beam spot. To form 7r° candidates, pairs of photon candidates with an invariant 
mass within 2.5a of the nominal 7r° mass are kinematically fitted with the mass constrained to the 
known ir° mass [1C]. 

The primary means of identification of B meson candidates is through their measured mass 



and energy. The beam-constrained mass of the candidate is defined as Mb = — |p| 2 , where 
p is the measured momentum of the candidate and Ef, is the beam energy. The resolution of Mb 
ranges from 2.5 to 3.5 MeV, where the larger resolutions correspond to decay modes with neutral 
pion(s). The second observable AE is defined as AE = E\ + E<i — E^, where E\ and Ei are the 
energies of the two final state mesons. The resolution of AE is mode dependent. For final states 
without a neutral pion, the AE resolution is about 20 MeV. For decay modes with one or two 
energetic neutral pions (B° — > n^p^, B° — > 7r°p° and B~ — > 7r°p~ etc), the AE resolution worsens 
by approximately a factor of 2 or 3 and becomes slightly asymmetric because of energy loss out of 
the back of the Csl crystals. We accept events with M B > 5.2 GeV and \AE\ < 100 to 300 MeV 
depending on the decay mode. 

The vector meson p, K* and u candidates are required to have masses within 200, 75 and 50 
MeV of their known masses [fLOf , respectively. In the simultaneous analysis of B° — > and 
tt°K* , the p° or K*° candidate is required to have mass between 0.3 GeV to 1.0 GeV under the 
7r + 7r~~ decay hypothesis so that both p° and K*° enter into the sample. Because of the polarization 
of the vector meson, the soft decay product from the vector meson may have momentum as low 
as 150 MeV. To reduce the large combinatoric background from soft 7r°s, only half of the helicity 
(H.) range, corresponding to a hard tt°, is selected when a p + or K* + decays to a tt + tt° or K + ir°. 
The helicity is defined as the cosine of the angle between one of the vector meson decay products 
in the vector meson rest frame and the direction of the vector meson momentum in the lab frame. 

The main background comes from continuum e + e~ — > qq, where q = u, d, s, c. This background 
typically exhibits a two-jet structure and can be reduced with event shape criteria. We calculate the 
angle 9s (0t) between the sphericity axis Q (thrust axis [l^]) of the candidate and the sphericity 
axis (thrust axis) of the rest of the event. The distribution of cos 0s(0t) should be flat for B mesons 
and strongly peaked at ±1.0 for continuum background. We require | cos^sl < 0.8 when there is 
a p or K* meson in the final state, and | cos 9t\ < 0.8 when there is aw meson in the final state. 
We also form a Fisher discriminant (J-) with event shape observables Q. 

We then perform unbinned maximum-likelihood fits where the likelihood of an event is pa- 
rameterized by the sum of probabilities for all relevant signal and background hypotheses, with 
relative weights determined by maximizing the likelihood function (£) The probability of 

a particular hypothesis is calculated as a product of the probability density functions (PDFs) for 
each of the input observables. The observables used in the fit are AE, Mb, T , Ti and the invari- 
ant mass of the resonance candidate. For final states with the same vector meson but different 
charged light mesons (pion or kaon), we also use the dE/dx measurement of the high-momentum 
track and fit for both modes simultaneously. Similarly, dE/dx measurements of the vector meson 
decay daughters are used in the simultaneous fit for B° — > tt°p° and 7r°i^*°. For each decay mode 
investigated, the signal PDFs are determined with fits to GEANT-based simulation |l3j samples. 
The parameters of the continuum background PDFs are determined with similar fits to simulated 
continuum samples as well as continuum data. Simulated continuum distributions are in excellent 
agreement with the data taken below the BB threshold. Correlations between observables used in 
the fits are investigated and their effect is found to be negligible. 

In all cases, the fit includes hypotheses for signal decay modes and the dominant continuum 



5 



background. Using the PDFs formed by the above observables, signal and continuum background 
can be well separated. For a few channels where the selected sample contains contributions from 
other B decays, we also include hypotheses for background from other B decay modes. These 
background decay modes can also be separated efficiently from the signal decay modes. We select 
a sample that contains both B — > vr^p , K~p° and some contamination from B~ — > tt~K*°. We 
then fit simultaneously for B~ — > vr^p , K~p° with and without a B~~ — > ir~K*° contribution. 
Similarly, we select a sample that contains both B~ — > ir~K*°, K~K*° with some contamination 
from B~ — > 7T _ /5°, K~p°. Then we perform a simultaneous fit for B — > it~~K* , K~K*° with 
or without the B~ — > 7r~p°, K~p° contributions. In both cases the fits with and without the 
background modes are consistent with each other. For each of the combinations i? — > 
B° — > 7r ± /3 =F , K+p*, and B~ — > 7T u, K u), contributions from other B decays are negligible and 
we select a common sample to fit for both modes. Finally individual samples are selected and fit 
for the B~ -» ir°p-, B~ -> tt°K*~, B° -> it uj and B° -» K° s uj searches. 

The contributions of b — > c and other i? decays are small in the selected samples of final 
states containing three tracks or two tracks and a 7r°, and their effects on the signal yields are 
negligible, except in the samples of B~ — > tt~p°, K~p° and B~ — > tt~K*°, K~K*°. Events from 
BT -» D°ir- where D° -> 7r T ,7r + 7r can enter into these samples and mimic our signal. We 
therefore impose a 30 MeV (~ 4<r) wide D° — > 7r + 7r" 7r T invariant mass veto in all the charged 
track pair combinations. We have also studied background from B~ — > K~rj', with r/' — > p°7 |l^ , [l0| . 
This background has exactly the same final state particles as -B - — > K~p° with an extra photon. 
Approximately 3% of this background can pass the selection for the £? - — ► vr~/9 , K~p° sample, 
therefore we include a component in the fit to describe this contribution. For B~ and 
5 _ — > tt°K*~ modes, due to the limited AE resolution for the final state with two neutral pions, 
the selected sample may contain background from other B processes such as B — > wax, pp. 

Table || shows the results of these measurements. The one standard deviation statistical error 
on the yield is determined by finding the ranges for which the quantity x 2 = — 21n£ changes by 
one unit. We observe significant yields for the decays B~ — > n~p , B° TC+p*, B- -nr w and 
B~ — >• 7r°/9 _ . To verify that the yields we observe in B meson decays to three-pion final states are 
indeed due to np decays, we repeat the standard fit allowing for an additional three-pion "non- 
resonant" contribution. The PDFs for this contribution are identical to the ones used for B — > irp 
signals except that we use PDFs that are constants in the p mass and TL. We find that this has 
no effect on the yield and the significance for B~ — > tt~ p° and B° — ► ir^p^ signals. Possible 
contributions from all other B processes, including higher mass pseudoscalar-vector decays, were 
also investigated for these channels and found to be negligible. However, the signal yield for 
B~ — > 7r p~ drops from 23.7^5' 4 with a significance of b.la to 8.0^y'g events with a significance of 
only la. We can not rule out the possibility that a significant fraction of the observed yield in ir°p~ 
comes from poorly measured processes such as non-resonant 7r _ 7r 7r°, irai and pp processes [10]. 
Therefore we calculate a conservative upper limit on the branching fraction assuming that the 
observed yield is due to B~ — > ir p~ decays only. 

Fig. 1 shows the likelihood contours from fits to B~ — > vr^p , K~p°, B° — > Tr^p^, K^p^ and 
B~ — > ir~u), K~lo. The resulting branching fractions are given in Table I. Fig. 2 shows the Mb and 
AE distributions after further requirements are made on event probability to reduce background. 
For the remaining processes in Table I we do not consider the signal yields to be significant (i.e. 
significance drops to less than 3 a after all the possible systematics are taken into account), and 
therefore set 90% C. L. upper limits for their branching fractions. Note that for the B~ — > K~to 
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TABLE I. Measurement results. Displayed are the decay mode, event yield from the fit, total 
efficiency including secondary branching fraction e, statistical significance (a), branching fraction 
from the fit Bfn (in units of 10~ 6 ), the measured branching fraction (B) or 90% confidence level 
upper limit (in units of 10~ 6 ) and theoretical prediction || (in units of 10~ 6 ). For the branching 
fraction measurement, the first error is statistical and the second systematic. We assume equal 
branching fractions for T(AS) — > B°B° and B + B~ . 



Decay 


mode 


Yield 


e(%) 


Signif. 




B or 90% B UL 


Theory 


B~ -> 


7T-p° 


29.8±|:| 


30 


5.4 


10.4l||±2.1 


10.4l|l ± 2- 1 


0.4 - 13.0 


B~ -> 


K-p° 




28 


3.7 


8A±f'° 4 ± 1.8 


< 17 


0.0 - 6.1 


B~ -> 




13.4i|;l 


18 


3.6 


7.6±1;^ ± 1.6 


< 16 


3.4 - 13.0 


B~ -> 


K-R*° 


o-o±g : g 


17 


0.0 


n n +i.3+o.6 
u - u -o.o-o.o 


< 5.3 


0.2 - 1.0 


B° — > 


■K±p* 


31.01H 


12 


5.6 


27.6±f:f ±4.2 


27.6l^; 4 ± 4.2 


12 - 93 


£?° — > 




1 fi 4+ 7 - 8 
iD -^-6.6 


11 


3.5 


16.01&2 ±2.8 


< 32 


0.0 - 12.0 


B° — > 


7T°p° 


k 4 +6.5 
°- 4 -4.8 


34 


1.2 


1.61^ ±0.8 


< 5.5 


0.0 - 2.5 


£?° — > 


7T°R*° 


o o +3 -° 

u - u -o.o 


25 


0.0 


n n+i.3+0.5 
u - u -o.o-o.o 


< 3.6 


0.7 - 6.1 


B" -> 


7T°p- 


23.7±f"t 


10 


5.1 


See text 


< 43 


3.0 - 27.0 


B~ -> 


TT°K*- 


2 6+ 4 - 2 

Z.D_2.6 


4 


1.0 


7.li^ 4 ± 1.0 


< 31 


0.5 - 24.0 


s- -> 


it~lj 


no c+8.2 
Zo.O_7 3 


26 


6.2 


11.3l|;jj ± 1.4 


11.3±1;| ± 1.4 


0.6 - 24.0 


£- -> 


K~U3 


7 q+6-6 
' - y -4.7 


26 


2.1 


3.2^ 4 ±0.8 


< 7.9 


0.2 - 14.0 


5°^ 


TT LU 


i c+3.5 


19 


0.6 


n q+1.9+1.0 

u -°-0.8-0.8 


< 5.5 


0.0 - 12.0 


5°^ 


K°LJ 


7 0+ 3 ' 8 

' - u -2.9 


7 


3.9 


10.01^ ± 1.4 


< 21 


0.0 - 17.0 



decay mode the additional CLEO II. V data and the re-analysis of CLEO II data no longer support 
its previously reported observation ||. However, the combined branching fraction B{B~ — * h~ui) 
= (I4.3I3 2 ± 2.0) x 10~ 6 (where h = K or ir) is still consistent with the previous result. 

Systematic errors are separated into two categories. The first consists of systematic errors in 
the PDFs, which are determined by varying the PDF parameters within their uncertainty. The 
second consists of systematic errors associated with event selection and efficiency factors. These 
are determined with studies of independent data samples. For branching fraction central values, 
the systematic error is the quadrature sum of the two components. For upper limits, the likelihood 
function is integrated to find the yield value that corresponds to 90% of the total area. The 
PDF systematic errors are taken into account in this procedure. The selection efficiency is then 
reduced by one standard deviation when calculating the final upper limit. As a goodness-of-fit 
check we compare — 21n£ at the minimum for our fits with expectations from fits to Monte Carlo 
experiments, and find them to be consistent in all cases. 

In summary, we have made first observation of the decays B~ — > vr _ p°, B° — > n^p^ and 
B~ — > ir~u). All of these AS = decay modes are expected to be dominated by hadronic b — > u 
transitions. We see no significant yields in any of the AS = 1 transitions. This is in contrast to 
the corresponding charmless hadronic B decays to two pseudo-scalar mesons (B — > PP) B — > Kir, 
7T7T, where AS = 1 transitions clearly dominate |l5]| . It indicates that gluonic penguin decays 
play less of a role in B — > PV decays than in B — > PP decays. This is consistent with theoretical 
predictions || that uses factorization which predicts destructive (constructive) interference between 
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penguin operators of opposite chirality for B — > X/? (B — > if7r), leading to a rather small (large) 
penguin contribution in these decays. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent luminosity 
and running conditions. This work was supported by the National Science Foundation, the U.S. 
Department of Energy, the Research Corporation, the Natural Sciences and Engineering Research 
Council of Canada, the A. P. Sloan Foundation, the Swiss National Science Foundation, the Texas 
Advanced Research Program, and the Alexander von Humboldt Stiftung. 
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FIG. 1. Likelihood contours at n standard deviations (a) of branching fractions for B — > it p°, 
K-p° (a), B° -» vripT, K±p* (b) and B~ -> ir~u, K~u (c). 
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FIG. 2. Projection plots in Mb and AE for B~ — > n p 
B~ — > tt~lj (e,f). The histograms show the data while the solid lines represent the overall fit 
to the data scaled to account for the extra requirement on event probability applied to make the 
projection. The dashed lines represent the continuum and the dotted lines on top of the contin- 
uum represent the other B components (B~ — > K~p°, B° — > K^p^ and B~ — > K~uj) in the 
simultaneous fits. 
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